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Abstract  -  A  useful  function  for  analyzing  burn-in,  developing  maintenance 
policies,  or  simply  modeling  lifetimes  of  equipment  is  the  mean  residual  life 
function.  Other  functions,  such  as  the  reliability  or  the  failure  rate 
functions,  are  of  course  important  also.  Discrete  data  arises  naturally  in 
various  ways:  from  discretizing  or  grouping  continuous  data,  devices  operate  by 
"cycles"  (e.g.,  a  copier's  "cycle"  is  a  copy,  its  "lifelength"  the  total  number 
of  copies),  etc.  This  paper  develops  a  general  approach  to  modeling  discrete 
bathtub  and  upside  down  bathtub  mean  residual  life  functions.  Because  the 
approach  allows  parametric  modeling  of  the  mean  residual  life,  maximum 
likelihood  estimation  of  models  can  be  done.  This  will  enable  estimation  of 
such  parametric  models  for  complete  discrete  data,  as  well  as  right  censored 
discrete  data.  A  simple,  perhaps  surprising,  example  is  presented  where  the 
mean  residual  life  increases,  then  decreases;  however,  the  hazard  rate  also 
increases,  drops  suddenly  at  one  cycle,  then  increases  again.  We  discuss  two 
reasonable  industrial  explanations  of  such  unusual  behavior. 


1 .  INTRODUCTION 


•>- 


Modeling  of  the  aging  process  of  components  and  equipment  can  be  performed 
in  various  ways.  One  helpful  tool  for  that  modeling,  as  well  as  for  analyzing 
burn-in  or  developing  maintenance  policies,  is  the  mean  residual  life  function. 
Of  course,  other  functions,  such  as  the  reliability  or  the  failure  rate 
functions,  are  important  also. 

An  introduction  to  MRL  (this  and  other  notations  are  defined  in  Section  2) 
and  its  many  applications  are  surveyed  in  Guess  and  Proschan  [4],  See  also 
Bhattachar jee  [2]  and  Gupta  [61.  Kuo  [9]  presented  a  review  in  his  Appendix  1 
on  burn-in  and  MRL.  Park  [10]  commented,  "...  the  time  at  which  a  bathtub 
failure-rate  is  a  minimum  does  not  maximize  the  mean  residual  life.  Contrary  to 
popular  belief,  the  mean  residual  life  in  the  constant  failure-rate  region  of  a 
bathtub  failure-rate  curve  is  not  constant."  Of  course,  this  and  other  MRL 
facts  have  important  consequences  for  developing  better  burn-in  policies. 

Discrete  data  arises  naturally  in  various  situations:  from  discretizing 
or  grouping  of  continuous  data,  devices  operate  by  "cycles,"  etc.  For  example, 
a  copier  and  its  key  components  have  as  their  "cycle"  producing  a  copy.  The 
lifelength  would  be  the  total  number  of  copies. 

Recently,  Ebrahimi  [3]  has  proposed  models  for  discrete  DMRL  and  IMRL.  He 
also  referenced  earlier  work  on  other  types  of  discrete  failure  models  for  life 
testing  data. 

Monotone  aging  models  are  very  useful  and  important  in  reliability 
applications.  E.g.,  a  Weibull  with  a  shape  parameter  greater  than  1  is  a  DMRL 
(also  an  IFR)  model  -  adverse  aging.  A  Weibull  with  a  shape  parameter  less  than 
1  (and  greater  than  0,  of  course)  is  an  IMRL  (also  a  DFR)  model  -  beneficial 


•\  ~W- 

•j  •_*  V  •j' 


.*vwv; 


aging.  Recall  IFR  implies  DMRL,  while  DFR  implies  IMRL.  On  the  other  hand, 
models  with  more  than  one  stage  of  aging  are  also  important  and  useful  in 
reliability  practice.  E.g.,  the  empirically  observed  bathtub  failure  rate  model 
has  aging  stages  of  infant  mortality,  useful  life,  and  wearout.  See  Barlow  s, 
Proschan  [1,  pp. 55-56]  for  an  example  on  a  commercial  airline  engine's  hot  gas 
generating  subsystem.  We  use  the  same  convention  as  Barlow  &  Proschan  [1,  p . 6 ] 
and  others  of  using  "increasing"  for  "nondecreasing"  (constant  is  allowed)  and 
"decreasing"  for  "nonincreasing." 

Note  that  bathtub  failure  rate  models  are  in  the  nonparametric  class  of 
DIFR.  The  related  class  to  DIFR  is  IDMRL.  The  lognormal,  used  for  repair 

times  as  well  as  lifetimes,  is  in  the  IDFR  class.  The  analog  to  I DFR  is  DIMRL. 
See  Guess,  Hollander  &  Proschan  [5]  for  details  on  the  IDMRL  and  DIMRL  classes. 

An  upside  down  bathtub  MRL  is  in  the  IDMRL  class  (human  lifelength  can  be 
modeled  well  by  this  class).  A  bathtub  MRL  is  in  the  DIMRL  class  (repair  time 
models  are  included). 

In  this  paper,  we  consider  discrete  versions  of  IDMRL  and  DIMRL.  A 
general  approach  to  creating  parametric  and  nonparametric  models  in  these 
classes  is  presented.  Illustrative  parametric  examples  are  given.  This  example 
includes  as  a  special  case  one  of  Ebrahimi's  [3]  model  by  allowing  n^  =  0. 

Because  the  approach  allows  parametric  modeling,  maximum  likelihood 
estimation  of  the  models  can  be  done.  This  estimation  would  be  possible,  not 
only  for  complete  discrete  life  data,  but  even  for  right  censored  discrete  data. 
Cf.  Hall  &  Wellner  (7). 

A  simple,  perhaps  surprising,  example  is  presented  of  IDMRL  but  not  DIFR. 
In  fact,  the  failure  rate  and  the  MRL  both  increase  over  the  set  n=0,l,...,  nQ-1 
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of  cycles.  At  the  turning  point,  nQ,  from  IMRL  to  DMRL  the  failure  rate  drops 
suddenly  and  increases  again.  See  Figures  1  and  2.  We  discuss  two  possible 
industrial  explanations  of  such  unusual  behavior.  Another  simple,  interesting 
example  shows  that  even  if  a  distribution  is  DIFR  and  I DMRL  the  turning  points 
in  each  can  differ.  See  Section  3. 


Notation 


I  DMRL 


DIMRL 


DIFR 


2.  NOTATION  &  NOMENCLATURE 


time  to  failure  of  a  system  (or  component) 


pmf  C  t i }  =  Pr{T=ti} 
Sf { t±}  =  Pr{T>t.} 
Sf{x}  =  Pr{T>x} 


discrete  failure  rate;  f./R. 

l  i 


complement  of  the  discrete  failure  rate;  s  1-X 

i  i 

mean  residual  life;  E(T-ti|T>ti) 
mean  residual  life 
decreasing  MRL 
increasing  MRL 

increasing  initially  then  decreasing  MRL 
decreasing  initially  then  increasing  MRL 
increasing  failure  rate 
decreasing  failure  rate 

decreasing  initially  then  increasing  failure  rate 
increasing  initially  then  decreasing  failure  rate 


Nomenclature 


I DMRL :  A  life  distribution  is  called  I DMRL  at  n„  if  m. <m.  ,  for 

0  i-  l+l 


vs 
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i=0,l,...,  n.-l  and  nu  >  mi+1  for  i=ng,  nQ+l,...  .  Note  that  the  MRL  is 
maximized  at  time  n^,  the  "turning  point." 

DIRML:  A  life  distribution  is  called  DIMRL  at  n q  if  nh  >  nu+^  for 
i=0,l, . . . ,Uq-1  and  nu  <  m^+1  for  i=ng,nQ+l, . . .  .  Note  that  the  MRL  is  minimized 
at  time  n^,  the  "turning  point." 

3.  A  GENERAL  APPROACH  FOR  CREATING  IDMRL  AND  DIMRL  MODELS 


Consider  again  the  definition  of  the  MRL  function 

m.  =  E[T-ti|T>ti]. 

This  can  be  rewritten  as 
00 

m i  =  J  R(x)dx/Ri 


(2.1) 


(2.2) 


Cf.,  for  example,  Kotz  &  Shanbhag  [8],  When  continuous  data  is  grouped  or  when 
it  is  only  possible  to  observe  failure  that  occurred  in  some  interval,  the  t^'s 
could  have  a  fractional  part  and/or  be  spaced  unevenly.  (E.g.,  during  night 
shifts  and  weekends,  the  inspection  times  for  failure  of  a  device  can  be  less 
frequent  and  sometimes  uneven.)  Allowing  for  this  and  assuming  the  t^'s  are 
ordered  (tQ  <  ^  <  ...),  we  can  express  (2.2)  as 


r  ”  r  j+1  1 

[U  =  E  J  R(x)dx  /TL 

lj=i  ;t.  J 

-  (.*  (tj+i  -  tj>  'h 


“j+i  y  ‘"'•j+i'j  '  i 

05 

'  C,fi  ">i  ’  h’Vi)  'h  l2-3) 

Recall  in  the  above  that  R(x)  is  left  continuous,  not  right  continuous. 

To  allow  for  fractional  values  and  uneven  spacings,  a  person  would  want  to 


use  (2.3).  On  the  other  hand,  for  the  case  t^+1  -  tj  =  1,  typical  of  devices 
(e.g.,  copiers)  with  cycles,  (2.3)  has  the  simple  form  of 


mi  "  .f.  Rj+1  7  Ri 


(2.4) 


A  simple  quick  check  of  (2.4)  is  to  take  T  degenerate  0.  (Note  that  Ro=l,  R^=0 
for  i>l ;  thus,  m0=Rx/Ro  +  0  +  0  +•••  =  0.  nu  =  0  for  i>l  by  convention). 

Because  we  have  been  motivated  primarily  by  devices  with  cycles  and 

previous  authors  have  focused  on  the  case  of  tj  =  j  for  j  =  0,1,2, we 

consider  (2.4)  for  the  rest  of  this  paper.  We  also  assume  nu  >  0  for 
i-0,1,2,...  in  all  of  the  following. 


The  main  tool  in  our  method  for  producing  IDMRL  and  DIRML  models  is  the 


following  theorem.  See  the  Appendix  for  its  proof. 


Theorem  1.  A  life  distribution  is  IDMRL  at  ng  if  and  only  if  there  exist 
a  sequence  { an } ,  an  >  0  for  all  n  such  that: 

an  *  an+l  for  H'0'1 . V1  |  (2  5) 

an  >  an+1  for  n=n0,nQ,+l,...  J 

and  Xn  =  an/(l+an+1)  for  n=0,l,2,...  .  (2.6) 

A  similar  result  for  a  DIMRL  at  ng  distribution  is  obtained  by  reversing 
the  inequalities  in  (2.5).  See  the  end  of  the  Appendix  for  further  comments. 


The  proof  in  the  Appendix  demonstrates  the  helpful  insight  that  an  s  mn. 
This  makes  building  models  even  easier.  In  creating  an  IDMRL  at  ng  model  we 
simply  need  to  generate  a  sequence  with  an  >  0  satisfying  (2.5)  and  (2.6).  Note 
that  condition  (2.5)  is  not  difficult  to  meet.  In  the  case  of  checking  whether 
a  known  distribution  is  IDMRL  at  ng,  then  (2.6)  is  also  straightforward. 
Condition  (2.6)  may  seem  rather  awkward,  however,  to  verify  in  the  situaton  of 
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creating  a  new  distributional  model.  Instead  of  verifying  (2.6)  directly,  it 
suffices  to  determine  that 
"  1 

E  —  =  ®  , 

n=0  an 

an+l  -  an  *  -1  '  n=0  > 1 / • •  • 

and,  of  course,  an  >  0  for  all  n.  When  a  sequence  is  increasing  after  some 
point,  the  above  inequality  is  always  met,  while  the  equality  may  or  may  not 
hold.  When  a  sequence  decreases  anywhere,  the  inequality  implies  the  decrease 
can  not  be  too  quickly.  When  a  sequence  is  decreasing  after  some  time,  the 
equality  will  always  hold.  (Aside:  it  turns  out  these  alternative  conditions 
are  necessary  ones  and  can  not  be  dropped.  They  will  always  be  implicit  in 
other  conditions.  Cf.  Bhattacher jee  [2]  and  Hall  &  Wellner  [7].) 

Although  the  increase  (decrease)  of  the  MRL  might  suggest  the  decrease 
(increase)  of  the  hazard  rate  of  the  life  distribution,  the  following  example 
shows  that  the  turning  point  for  IDMRL  does  not  necessarily  agree  with  the 
turning  point  from  decreasing  hazard  rate  to  increasing  hazard  rate.  Compare 
also  Example  2. 

Example  1.  Let  aQ  =  1,  a-^  =  3/2,  a2  =  4/3,  a^  =  1  and  an  «=  1/3  for  n  >  4. 
Then  it  is  easy  to  verify  that  the  turning  point  from  DFR  to  IFR  occurs  at  n=3. 
However,  according  to  Theorem  1,  the  turning  point  for  IDMRL  occurs  at  Dq=1. 

Ebrahimi  [3]  considered  the  following  illustrative  discrete  DMRL  model 


an  a  c/(«n+l),  n=0,l,2,...,  (0<c<l,  a>0). 


(2.7) 


It  can  be  shown  directly  that  an  =  mn  <  1  for  n=0,l,2,...  for  model  (2.7).  Note 
that  mn  <  1  is  a  severe  restriction,  of  course.  A  very  slight,  natural 
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modif ication  of  the  model,  however,  is  possible.  Let  T  be  a  random  variable 

with  distribution  formed  from  the  (a  }  of  (2.7).  The  sequence  {a  )  that 

n  n 

f 

generates  the  distribution  of  T  =  9T  for  0>O  is 
2 

an  =  0  c/(an+0),  n=0,l,2,...,  (0<c<l,  a>O,0>O).  (2.8) 

Note  that  E(T  )  =  0*E(T)  and  that  for  0=1  (2.8)  reduces  to  (2.7).  Also  recall 

m  s  a  . 
n  n 

To  get  a  simple  discrete  IDMRL  model,  we  let 
an  =  ©2c/(oc|n-n0|  +  0) 

for  n=0,l,2,...,  (0<c<l,  a>0,  0>O,r.Q>O).  An  easy  check  of  Theorem  1  shows  the 
above  valid.  This  model  includes  both  (2.7)  and  (2.8)  as  special  cases.  (The 
case  nQ=0  is  allowed.) 

It  is  possible  to  simplify  this  IDMRL  an  by  letting  c  =  1;  thus, 
an  =  02  /  (a  |  n-nQ|  +  0). 

Note  that  lim  an  =  0  is  a  natural  condition  on  the  MRL.  Recall  the  failure 
n-*® 

rate  of  the  Gamma  life  distribution,  however,  converges  to  a  constant  greater 

than  0  as  time  t  gets  larger.  I.e.,  asymptotically  in  time  it  behaves  like  an 

exponential.  If  a  situation  suggests  that  lim  a  =  a  >  0  is  needed  on  the  MRL, 

n-x» 


another  possible  modification  is 

an  =  02/(a|n-nQ|  +  0)  +  y,  n=0,l,2,...  ,  (2.9) 

with  the  new  parameter  y  >  0.  If  a  reliability  analyst  wanted  the  essential 
flexibility  of  (2.9),  yet  a  simplified  version  then  just  set  a  =  1. 

Example  2.  Let  a  =  1,  n^  =  40,  0  =  62.711168,  and  y  =  61.711168  in  (2.9).  To 
get  a  mean  ( m g)  of  100,  we  used  this  0  and  y.  The  MRL  function,  nu,  is  graphed 
in  Figure  1.  Note  that  the  graph  is  as  expected  IDMRL  at  ng  =  40.  Figure  2 


shows  the  graph  of  the  failure  rate  function.  Note  that  it  is  not,  however, 
DIFR.  In  fact,  both  the  failure  rate  and  the  MRL  strictly  increase  for  n=0, 
1,...,  39=Oq-1.  (Recall  IFR  everywhere  implies  DMRL  everywhere . )  At  the 
turning  point  (ng  =  40)  from  I MRL  to  DMRL  the  failure  rate  drops  suddenly  to  its 
minimum  and  then  increases  again. 

One  possible  explanation  for  the  unusual  behavior  is  that  at  the  turning 
point  an  overhaul  is  performed  with  new,  state  of  the  art  parts.  The  higher 
initial  IFR  is  replaced  by  a  lower  IFR  immediately  upon  completion  of  the 
overhaul.  The  closer  in  time  the  device  can  get  to  the  overhaul  time  and  the 
significantly  lower  failure  rate,  the  larger  the  MRL  grows.  This  explains  the 
initial  I MRL  part,  even  though  the  failure  rate  is  increasing.  In  the  tail 
(n>ng  =  40)  we  have  IFR,  thus  DMRL. 

Another  explanation  is  to  consider  a  stress  screening  program  of  parts  or 
devices  with  IFR.  in  the  initial  time  period  the  stressful  screen  results  in  a 
much  higher  than  normal  IFR.  This  is  followed  upon  completion  by  a  sudden  drop 
to  a  significantly  lower  IFR.  The  closer  we  get  to  the  time  of  exiting  the 
screen,  the  higher  the  MRL  gets,  i.e.,  we  have  I MRL  initially.  After  exiting  to 
the  lower  IFR  we  have  DMRL. 

Figures  1  and  2  also  show  the  importance  of  graphing  both  the  MRL  and  the 
failure  rate  to  understand  data  better.  It  will  not  necessarily  be  obvious 
looking  at  one  graph  what  the  other  does. 

4.  CONCLUSIONS 

For  cost  analysis,  burn-in,  and/or  modeling  discrete  life  data,  the  MRL  is 
a  useful  tool.  We  have  provided  a  general  approach  for  creating  I DMRL  (e.g., 
useful  for  lifetimes)  and  DIMRL  (e.g.,  useful  for  repair  times)  models. 
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Given  a  parametric  version  of  one  of  these  models,  it  is  possible  using 
(2.6)  to  derive  the  likelihood  for  the  data.  See  Ebrahimi  [3]  and  his 
references.  Cf.  Hall  and  Wellner  (7j.  With  the  likelihood  it  becomes  possible 
to  do  maximum  likelihood  estimation  on  complete  discrete  life  data  and  even  on 
right  (and/or  left)  censored  data. 

An  advantage  to  parametrizing  directly  the  MRL  as  presented  in  this  paper 
is  that  even  given  the  failure  rate  or  the  probability  mass  function, 
calculation  of  the  MRL  could  require  doing  an  infinite  summation  numerically. 


and 


f. 

1 


i-1  _  i-1 

\  n  X  =  [l-a./<l+a .  .)]  n  [a/d+a..)] 
1  j=0  3  1  1+1  j=0  3  3+1 

for  i  =  1,2, .. . 


i-1  _  i-1 

R.  =  n  X.  =  n  [ a  ./{ 1+a  •  , ) ] 
j=0  3  j=0  3  3+1 

for  i=l,2, . . . 


where  {an}  is  a  sequence  defined  in  Theorem  1. 

To  prove  Theorem  1,  we  need  the  following  lemma  from  Ebrahimi  [3]. 


Lemma  1.  Let  {an}  and  {b^}  be  two  converging  positive  sequences  and  let  c  >  0. 

If  ar/^c+an+l'>  3  kr/^c+l:)n+l^  f°r  a3-3-  n»  then  an  s  bn  for  all  n. 

Proof.  The  proof  of  the  lemma  involves  essentially  case  I:  a  s  b  t  0  and  case 

II:  a  =  b  =  0,  where  lim  a„  s  a  and  lim  b  e  b. 

n-*»  n  n-»®  n 


Proof  of  Theorem  1.  First  assume  that  the  life  distribution  is  IDMRL  at  n^. 
Recall  we  have  assumed  mn  >  0  for  all  n.  Note  that 


mn  ‘  A  V^n  ■  ^ 

-  .  *  l;1  <V/V 

i=n+l  j=n  J  J 
<*>  i-1 

=  z  n  x. 

i=n+l  j«*n  3 

®  i-1 

-  X  +  X  Z  n  X. 

n  n  .  _  . 

i=n+2  ]=n+l  J 


(A.l) 


-  Xn  (1  +  mn+^)  holds  for  any  MRL,  mn. 


Thus,  we  express 


Xn  .  mn/(l  +  mntl)  . 


(A. 2) 


By  definition  of  IDMRL  at  nn  the  sequence  {a  }  defined  by  a  =  m  satisfies 

v  n  n  n 

(2.5).  From  (A.2),  {a  1  also  fulfills  (2.6).  Note  that  a  =  m  >  0  for  all  n. 

n  n  n 

This  completes  the  first  implication  of  Theorem  1. 

To  prove  the  implication  in  the  other  direction,  we  assume  the  existence 

of  a  sequence  {a  1  satisfying  (2.5),  (2.6)  and  a  >  0  for  all  n.  By  a 
n  n  n 

decreasing  for  n  >  nQ  and  an  >  0,  lim  an  =  a  >  0  exists. 

n->® 

Now  consider 


oo  i_i 

l  II  X. 
i=n+l  j=n 


(A.  3) 


®  i-l+n 
h  r  n  X, 
i=l  j=n  J 


which  follows  directly  from  the  representation  of  X^  in  (2.6),  the  definition 
of  mn,  ( A. 1 ) ,  and  readjusting  the  summation  index.  Note  that 


lim  Xn  s  lim  a  /(1+a  .)  =  a/(l+a). 

n-*=  n-*” 


(A.  4 ) 


By  (A. 4)  and  the  ratio  test,  mn  converges,  and  hence,  is  well  defined. 


Consider 


lim  mn/(l+mn+1)  =  lim  t^n(1+mn+1)  ]/U+mn+1) 


=  lim  Xn  =  a/(l+a)  . 
n-*» 


(A.5) 


Using  limiting  operations,  we  have 
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®  (1-1+n 

lim  m  =  I  lim  n  X . 

n  .  .  .  i 

n-x°  i=l  n-*=°  ^  ]=n  J 


By  \n  =  an/(l+an+1)  =  n^/(l+mn+1)  and  the  results  above,  we  apply  Lemma  1  to  see 
an  =  mn  for  all  n.  By  property  (2.5)  of  (an) ,  the  sequence  {mn)  satisfies  the 
condition  for  IDMRL  at  ng,  which  we  wanted  to  show. 

For  the  DIMRL  at  ng  version  of  Theorem  1  the  proof  assuming  first  DIMRL 
(nQ)  follows  similarly.  To  show  the  implication  in  the  other  direction, 
however,  requires  the  additional  condition  that  the  sequence  {an}  is  bounded  for 
these  proof  techniques.  Watson  &  Wells  [12]  considered  lognormal's  unbounded 
MRL.  We  allow  for  unbounded  { an)  and  MRL  by  our  comments  and  conditions  on  page 
7.  Cf.  also  Bhattacharjee  (2]  and  Hall  &  Wellner  [7). 
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